function [post act] = gmmpostRobustVerbeek(mix, data)


K = mix.ncentres;
N = size(data.Y, 1);

[act logact] = gmmactivRobustVerbeek(mix, data);
log_priors = log(mix.priors);

logpost = repmat(log_priors, N, 1) + logact;
max_logpost = max(logpost, [], 2);
logpost = logpost - repmat(max_logpost, 1, size(logpost, 2));

post = exp(logpost);
post = post./(sum(post, 2)*ones(1, K));
